Antiferromagnetism in the magnetoelectric effect single crystal LiMnP04 
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Elastic and inelastic neutron scattering studies reveal details of the antiferromagnetic tansition 
and intriguing spin-dynamics in the magneto-electric effect single crystal LiMnP04. The elastic 
scattering studies confirm the system is antiferromagnetic (AFM) below rjv=33.75 K with local 
magnetic moments (Mn^^; S — 5/2) that are aligned along the crystallographic a-axis. The spin- 
wave dispersion curves propagating along the three principal axes, determined by inelastic scattering, 
are adequately modeled in the linear spin-wave framework assuming a spin-Hamiltonian that is 
parameterized by inter- and in-plane nearest- and next-nearest-neighbor interactions, and by easy- 
plane anisotropy. The temperature dependence of the spin dynamics makes this an excellent model 
many-body spin system to address the question of the relationship between spin-wave excitations 
and the order parameter. 

PACS numbers; 75.25.-|-z, 75.30.Ds, 75.50.Ee 
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INTRODUCTION 



The recent discoveries of colossal magnetoelectric ef- 
fects (ME) in rare-earth-manganites (RMnOa P, 0) and 
manganese-oxides (R2Mn05 [3|]) triggered a revival in- 
terest in the so-called insulating mulitiferroic materials 
that exhibit ferroelectricity in coexistence with ferro- 
magnetism or antiferromagnetism (FM or AFM) 0, @] ■ 
Systematic studies of the coupling between the electric 
and magnetic fields in crystals date back to the early 
1960s with the discovery of the first ME compound 
CraOa UM- Early on, the isostructural transition- metal 
lithium-orthophosphates Li7kfP04 (M = Mn, Fe, Co, Ni) 
were identified as ME systems [IJs 
the subjects of numerous studies [ll 



1C| and have been 
12, HI. Like other 



members of the lithium-orthophosphates, LiMnP04 is an 
antiferromagnetic insulator with Prima symmetry group 
14, [l^. In this structure, each Mn^+ ion occupies 
the center of a slightly distorted MnOg octahedron that 
shares oxygen anions with a tetrahedral PO4 forming 
a closely packed oxygen framework. The Mn^+ ions 
{S — 5/2) form buckled layers that are stacked along 
the [100] crystallographic axis, as shown in Fig. [Ija). 
The nearest neighbors (NN) in the b-c plane are coupled 
magnetically by a relatively strong exchange interaction 
Ji through an Mn-O-Mn oxygen-bond, whereas the in- 
plane next-NN (NNN) are coupled via Mn-O-O-Mn (J2) 
|16l. Il7| (see Fig. [Ijb) for the definitions of the exchange 
couplings). The interlayer magnetic coupling is mediated 
through phosphates by higher order superexchange via 
Mn-O-P-O-Mn, which was found to be relatively large in 
similar frameworks [isl. 



Neutron diffraction of polycrystalline samples [l9|, |20|, 
21 1 and single crystal NMR [l^| measurements showed 
that all LiAfP04 share the same collinear (up-down) 
AFM ground state with spin orientation along a, b, b 
and c crystallographic directions for LiMnP04, LiFeP04, 
LiCoP04 and LiNiP04, respectively. However, recent 
single crystal neutron diffraction studies of LiCoP04, 
LiFeP04, and LiNiP04, [H [H [H [H show the mo- 
ments in the ground state are slightly tilted away from 
principal crystallographic directions, indicating the mag- 
netic symmetries for these systems are lower than those 
determined from polycrystalline measurements, giving 
rise to spontaneously induced weak ferromagnetism. 
Weak ferromagnetism (WFM) in magnetic susceptibility 
measurements has also been reported for LiNiP04 |26[ 
and LiMnP04 27| below Tjy. Indeed, domain structures 
observed by second-harmonic-generation (SHG) experi- 
ments in LiCoP04 were interpreted as ferrotoroidic do- 
mains [2^ facilitated by the lower magnetic symmetry 
obtained in neutron scattering experiments [22*1. Based 
on the detailed spin configuration observed in LiNiP04, 
Jensen and co-workers have been able to model the tem- 
|)CTature dependence of the ME coefficients of this system 

Here, we report elastic and inelastic neutron scatter- 
ing studies of a single crystal LiMnP04, to determine 
the nature of the AFM transition and the spin dynamics 
in this system. Recent susceptibility measurements in- 
dicated WFM in this system [27[ implying spin-canting 
that may be detected in neutron diffraction measure- 
ments. There is also some inconsistency in the literature 
with regard to the transition temperature; T/v =34.85 
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20| and 42 K were reported [27[ for polycrystalline 
samples. The spin dynamics of the LiFeP04, LiCoP04 
and LiNiP04 were measured and modeled in the linear 
spin-wave framework only recently [22 



24, \m, to sue- 
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cessfuUy yield the exchange coupHngs and the single-ion 
anisotropy parameters in these systems. Determining 
and analyzing the spin dynamics of LiMnP04 is an im- 
portant step towards developing a universal understand- 
ing of the magnetic properties of this isostructural group 
of compounds. 





FIG. 1: (color online) (a) Atomic structure of LiMnP04. 
The Mn2+ ions form buckled layers stacked perpendicular 
to the [100] crystallographic direction. The ground state 
of LiMnP04 is coUinear antiferromagnetic with average mo- 
ments along the a-axis. (b) Spin arrangement of the two Mn^"*" 
layers. The in-plane nearest and next-nearest neighbor inter- 
actions Ji, J2, J3 and inter-plane nearest and next- nearest 
neighbor interactions J4, J5 are labeled. 



EXPERIMENTAL DETAILS 



(99.999% Aldrich) [30|. Powder, for the XRD, was pre- 
pared by crushing typical isolated single crystals. The 
composition and structure were confirmed by carrying 
out Rietveld analysis of X-ray powder diffraction (XRD) 
data, using the GSAS software package [3l|. No extra 
peaks from impurities were detected in the XRD pat- 
tern. The lattice parameters yielded from the refinement 
at room temperature (a = 10.524 A, 6 = 6.095 A, and 
c — 4.75 A) are in good agreement with the values re- 

[13,133]. 



ported in the literatures [2C 



Neutron scattering measurements were carried out on 
the BT7 and BT9 thermal triple axis spectrometer at the 
National Institute of Standards and Technology (NIST) 
Center for Neutron Research (NCNR) . A monochromatic 
neutron beam of wavelength A = 2.36 A (14.7 meV, 
/to = 2tt/X = 2.66A-1) was selected by a vertical fo- 
cusing monochromator system, using the (0 2) Bragg 
reflection of highly oriented pyrolytic graphite (HOPG) 
crystals. HOPG crystals were also used as analyzer for 
both the elastic and the inelastic studies. The high reso- 
lution inelastic scattering measurements were conducted 
on the on the cold neutron Spin Polarized Inelastic Neu- 
tron Spectrometer (SPINS) at the NCNR. 



RESULTS AND DISCUSSION 
Elastic neutron scattering 

The LiMnP04 crystal was oriented with its a-b plane, 
(and subsequently rotated with its b-c plane) to coin- 
cide with the scattering-plane of the spectrometer. The 
elastic measurements, with the strongest magnetic reflec- 
tion (010) peak, conflrmed that the magnetic structure of 
LiMnP04 is AFM with spin orientation along the a-axis. 
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A LiMnP04 single crystal (0.41 gram, pink in color) 
was grown by the standard flux growth technique (LiCl 
was used as the flux) from a stoichiometric mixture 
of high purity MnCla (99.999% Aldrich) and Li3P04 



FIG. 2: (Color online) Temperature dependence of the square 
root of the integrated intensity of the (010) magnetic peak. 
The transition temperature obtained from the fit is Tn = 
33.85 ± 0.1 K and the critical exponent /? = 0.126 ± 0.017. 
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The temperature dependence of the magnetic (010) re- 
flection was used to extract the behavior of the order 
parameter in the temperature range from 5 to 50 K. Fig- 
ure m shows the square root of the integrated intensity, 
representing the staggered magnetization, (i.e., order pa- 
rameter) as a function of temperature for the (010) peak. 
The order parameter was fit to a power law function near 
the transition temperature: 

V7 cx Af^ = M^tf^ (1) 

where Mq is the sublattice magnetization at T = K, 
t — {1-T/Tm) is the reduced temperature, and (3 is the 
critical exponent. The obtained transition temperatures 
from the fit is TV = 33.85 ± 0.1 K and the critical ex- 
ponent for the temperature dependent magnetization /3 
is 0.126 ± 0.017 by using the main (010) magnetic peak. 
This is very close to the theoretical value of the critical 
exponent of the 2D Ising system (3 =0.125 3J], consis- 
tent with the layered nature of the magnetic system, as 
also demonstrated by the weak interlayer coupling ob- 
tained from the analysis of the spin-waves discussed be- 
low. The transition temperature is found to be very close 
to the value, 34.85 ± 0.1 K measured by Mays [l^| using 
nuclear magnetic resonances performed on a single crys- 
tal of LiMnP04, whereas susceptibility measurements of 
powder samples yield Tn = 42 K [271]. 

Unlike LiCoP04 and LiFeP04, strong critical scatter- 
ing above the AFM transition is observed in LiMnP04 
and persists to almost twice Tjv (traced to tempera- 
tures as high as 70 K) before the spins become uncor- 
related. These correlations were already evident in Fig- 
ure m Figure [3] shows longitudinal and transversal scans 
at the (010) magnetic peak above the transition with 
energy transfer AiJ = 0. The peaks are much broader 
than the spectrometer's resolution indicating some type 
of short range correlations. This is reminiscent of the 



behavior in LiNiP04 where this critical scattering [21 1 
was later found to be associated with an incommensu- 
rate (IC) short- and long-range magnetic order above 
Tn [35|- Below ^33 K all peaks are practically resolu- 
tion limited Gaussian shaped. Above T/v, the peaks were 
fit to a Lorentzian line shape l/{q^ -t- k^), where q = h 
or fc, and k is inversely proportional to the coherence 
length ^ = 27r/K, convoluted with an instrumental gaus- 
sian shaped resolution function. The calculated coher- 
ence lengths along the a- and 6-axis, as a function of tem- 
perature are shown in Figure [3] (c) . Below the transition 
temperature, the in-plane coherence lengths (along the b- 
axis) is significantly longer than that between the planes 
(along the a-axis), consistent with the quasi-2D nature of 
LiMnP04. The correlation lengths were remeasured on 
BT9 using the integrated energy (two-axis mode) method 
and yielded similar results. Attempts to correlate the in- 
plane coherence length in the critical regime with the 
2D-Ising model failed, but the Kosterlitz-Thouless (KT) 
2D XY-model, C(T) = Ae^/C^"'^'^)" , seems to fit our data 
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FIG. 3: (Color online) (a) Longitudinal and (b) transver- 
sal, (i.e., inter-plane) critical magnetic scattering scans at 
the (010) magnetic reflection above Tn- (c) Correlation 
lengths obtained after deconvoluting the spectrometer's res- 
olution function (in-plane and inter-plane as indicated) ver- 
sus temperature. Solid line is calculated assuming Kosterlitz- 
Thouless theory. The error bars in this paper are statistical 
in origin and represent one standard deviation, (r.l.u. stands 
for reciprocal lattice unit) 



weh [361 (sofid fine. Fig. O with = 33.6 ± 0.008 and 
V = 0.51 ± 0.1. This strong critical scattering above the 
transition with KT characteristics may therefore indicate 
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spin-dimensionality crossover from the 2D Heisenberg to 
the 2D XY-model. 



Inelastic and quasielastic neutron scattering 

Spin waves along the three principal reciprocal lattice 
directions {q,l, 0), (0, 1+q, 0) and (0, 1, q) were measured 
in energy loss mode at T = 5 K. Examples of the excita- 
tions measured on BT7 at q — 0.2 are shown in Figure 
m A single excitation was observed at each q along the 
three directions on BT7 which has an energy resolution 
of ~1 meV. At T w 50 K, no similar peaks are observed 
confirming the magnetic origin of the excitations. The 
inelastic signals at various constant wave-vectors q were 
fit to Gaussian shaped functions (solid line in Fig. [4|), 
and the set of energies at maximum intensity were used 
to construct the spin-wave dispersion curves shown in 
Fig. [6l It is shown that the spin- waves propagating in the 
plane along the (001) and (010) directions have higher en- 
ergy than the spin- waves propagating along (100) at the 
same q values. Qualitatively, this behavior reflects the 
anisotropy in the strength of the exchange couplings in 
the system; as expected, the in-plane exchange couplings 
are much stronger than those between planes. Using 
the cold neutron triple axis SPINS spectrometer, an en- 
ergy gap Eg = 0.48 meV was observed around the (010) 
zone center, which is much smaller than the 2 meV p9| . 
5.86 meV and 4.7meV ^ observed in LiNiP04, 
LiFeP04 and LiCoP04, respectively. With the high en- 
ergy resolution of SPINS, which is around 0.1 meV (using 
3.7 meV final energy), two energy excitation peaks were 
identified at the zone center, as shown in Figure [51 
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FIG. 4: (color online) Examples of constant-Q energy scans 
measured on BT7 at 5 K, at wave- vectors q =0.2 along {q, 
1, 0), (0, l+q,0) and (0,1, g) reciprocal directions. A single 
energy excitation is present in every direction with the typical 
energy resolution of BT7 around 1 meV. 

To model the spin-wave dispersions, we use a spin 



Hamiltonian based on the ground state spin structure 
of LiMnP04 as shown in Fig. [1] as follows 

Ji to J5 are the spin couphng constants (see Fig. [TJ, 
and Dj are the single ion anisotropics. Since the exci- 
tation spectrum is insensitive to an overall shift of the 
ground state energy we can define = for simplic- 
ity. The x,y, and z coordinates are defined along the 
C-, b- and a- axis, respectively, to align the spin direc- 
tion in the ground state with the quantum z-axis in Eq. 

The magnon dispersion curves derived from Eq. ^ 
by linear spin- wave theory is given in Refs [23j, |29|, l37j . 
In the model, the calculated spin waves have two non- 
degenerate branches (denoted by the ± sign in Eq. [3] ) 
as a result of the different anisotropics along the x and y 
directions. 

huj = -(B± C)2 (3) 

where, 

A EE 45(Ji + J5) - 2^[J3(1 - cos(q • rs)) + J2(l - 
cos(q • rg)) + J4(2 - cos(q • rr) - cos(q • rs))] 
+ {S~l/2){D, + Dy), 

B={S-l/2)iD,-Dy), 

C= 2Ji5'[cos(q • ri) + cos(q • r2)] 

-|-2J5S'[cos(q • rs) -I- cos(q • r4)], 

and Ti denotes a vector to a NN and NNN, 

ri = (c/2, 6/2,0) r2 = (-c/2, 6/2,0) 

rg = (0, 6/2, a/2) r4 = (0, -6/2, a/2) 

r5 = (0,6,0) r6 = (c,0,0) 

r7 = (c/2,0,a/2) r^ = (-c/2, 0, a/2). 

The spin-wave dispersion curves along the three direc- 
tions in Fig. [5] were simultaneously fit to Eq. ([3]), using 
the "— " sign, yielding the following values: Ji = 0.48 ± 
0.05 meV, J2 = 0.2 ± 0.038 meV, J3 = 0.076 ± 0.004 
meV, J4 = 0.036 ± 0.002 meV, J5 = 0.062± 0.003 meV, 
^ 0.0069 ± 0.001 meV and Dy = 0.0089 ± 0.001 meV. 
In the equation, S* = 5/2 for Mn^+. As expected, the 
in-plane NN exchange coupling Ji, is the strongest, com- 
pared to the in-plane NNNs J2 and J3 . The sign of both 
J2 and J3 indicates the NNN interactions compete with 
the simple AFM ordering dictated by Ji. For weakly 
coupled layers, it has been predicted theoretically that 
an incommensurate (IC) majgnetic structure should be 
realized when J2/J1 > 0.5[3^ Thus, unlike in LiNiP04 
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where the J2/J1 « 0.6 and an IC has been observed, 
the ratio for LiMnP04 (~ 0.4) seems to be too small to 
induce any IC phase transition [2§|. The spin couplings 
between the inter-plane nearest-neighbors ( J4 and J5) are 
relatively weak at about 12% of Ji consistent with the 
quasi-2D behavior of this system. The values of the single 
ion Dx = 0.0055 and Dy = 0.0071j_are much smaller than 



24| 



those of LiNiP04 |29|, LiFeP04 [23l, and LiCoP04 
indicating that the ground state with magnetic moments 
along the a-axis is not very stable, and the moments are 
prone to a spin- flop transition in relatively weak magnetic 
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FIG. 5: (color online) The energy excitation at the zone center 
measured on SPINS cold neutron triple axis which has an 
energy resolution of ~ 0.1 meV. Two excitations are clearly 
identified at the zone center. 

The second spin wave dispersion branches, given by 
gjgjj jj^ gq are calculated using the Js and Ds 
obtained from the fits listed above. The two branches al- 
most overlap one another for the dispersions along all the 
three principal reciprocal directions, and are only sepa- 
rate by 0.1 meV at the zone center. The spin wave dis- 
persion along (g, 1, 0) direction, where the model predicts 
the largest separation between the two branches, was re- 
measured with the high energy resolution on SPINS. Fig- 
ure [H] (b) shows an enlargement graph of Figure [S] (a), 
with fairly good agreement with the model calculations. 

The energy gaps at the zone center for the two branches 
are 



A£; = 2S^AD4Ji + J5) + D^Dy, 
for (B - C) in Eq. © and 



AE = 2SJADy{Ji + J5) + D^D 



1" 



(4) 



(5) 



for (B -I- C). J5 represents the inter-plane NN coupling. 
From the equations, we notice that the energy gap not 
only depends on the single-ion anisotropy terms, but also 
on the two nearest-neighbor antiparallel exchange inter- 
actions. 
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FIG. 6: (color online) (a)Spin-wave dispersion curves along 
the a* , b* , and c* reciprocal space directions measured at 5 
K. Solid lines are best-fit calculations obtained from linear 
spin- wave theory using Eq. (b) Zoomed plot of (a) near 
the zone center. The predicted second spin wave dispersion 
branches are shown as dashed lines. 



0.8 

> 
CD 

^ 0.6 

Q. 
03 

O 0.4 
>. 

05 0-2 
c 
LU 

0.0 

















\ i 



10 15 20 25 30 35 

Temperature (K) 



FIG. 7: (Color online) The energy gap as a function of tem- 
perature measured on the SPINS cold neutron triple axis. 



The temperature dependent energy gap up to the tran- 
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sition temperature was measured at the cold neutron 
triple axis spectrometer SPINS, and the results are shown 
in Figure[7] (energy gaps at various temperatures were de- 
termined from gaussian fits to constant-Q energy scans 
such as the one shown in Fig. [5] at T = 5 K). The energy 
gap monotonically decreases with increasing temperature 
and approaches zero at the transition temperature. The 
temperature dependence of the gap to a first approxi- 
mation is proportional to the st agg ered magnetization 
which is temperature dependent [42] ■ However, it may 
deviate in the critical regime due to the different tem- 
perature dependencies of the coupling constants and the 
single ion anisotropy. In antiferromagnets, the exchange 
constants Jg usually decrease much faster than the s ing le 
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ion anisotropy near the transition temperature (4l|, |42 1 



Quasi-elastic scattering (QENS) around (010) at differ- 
ent temperatures was measured on BT9 using the inte- 
grated energy (two-axis) mode, and the results are shown 
in Figure [5] At temperatures right below the transition, 
the (010) peak consists of a resolution limited Gaussian 
shaped magnetic Bragg peak superimposed on a broad 
Lorentzian shaped diffuse peak. Whereas the diffuse scat- 
tering becomes stronger with the increase of temperature 
(up to the transition), the elastic magnetic Bragg peak 
becomes weaker. The QENS intensity at each tempera- 
ture was integrated over the the K range shown in Figure 
[8](a) excluding the region from 0.98 to 1.02 (r.l.u) which 
is dominated by elastic scattering. Figure [8] (b) shows 
the QENS as a function of temperature, which exhibits 
a sharp peak at the transition (T^v = 33.75 K) with a 
tail that extends to about 1.5Tn. This indicates that the 
short range correlations observed in the elastic scattering 
are primarily due to (dynamics) spin- fluctuations. 

In summary, we determined the critical behavior near 
the AFM magnetic phase transition of LiMnP04 {Tn = 
33.85 K). The strong critical scattering around the (010) 
magnetic peak and the in-plane, inter-plane coherence 
lengths indicate that the system is a quasi-2D system 
with very weak easy axis single ion anisotropy. Analysis 
of the spin- wave dispersions along the three principal axis 
directions show that the in-plane couplings are dominant 
compared to the inter-plane couplings. These in-plane 
competing interactions between in-plane NN and NNN- 
spins in LiMnP04 seem to be too weak to lead to more 
complicated, incommensurate magnetic structures. This 
is in contrast to the observation of incommensurate mag- 
netic phases in LiNiP04 [s^. 

We thank Sung Chang for the support on the SPINS 
cold neutron triple axis. Ames Laboratory is supported 
by the U.S. Department of Energy, Basic Energy Sci- 
ences, Office of Science, under Contract No. DE-AC02- 
07CH11358. SPINS is supported in part by the US Na- 
tional Science Foundation through DMR-0454672. 
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FIG. 8: (Color online) (a) Examples of the quasielastic scat- 
tering scans around (010) magnetic peak at temperatures be- 
low and above the transition temperature Tjv, which were 
measured on BT9 using the integrated energy method, (b) 
The temperature dependent integrated intensity from the 
quasielastic scattering excluding Q = 0.98 to 1.02 as indicated 
by the box in (a). 
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